Spontaneous Quantum Hall States in Chirally-stacked Few-Layer Graphene Systems 
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Chirally stacked iV-layer graphene systems with N > 2 exhibit a variety of distinct broken sym- 
metry states in which charge density contributions from different spins and valleys are spontaneously 
transferred between layers. We explain how these states are distinguished by their charge, spin, and 
valley Hall conductivities, by their orbital magnetizations, and by their edge state properties. We 
argue that valley Hall states have [N/2] edge channels per spin- valley. 
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PACS numbers: 73.43.-f, 75.76.+j, 73.21.-b, 71.10.-w, 75.85. +t 



Introduction — In the early 1980s, following the discovery 
of the quantum Hall effect (QHE)[l|, it was recognized [2| 
that electronic states can be characterized by topological 
indices, in particular by the integer valued Chern num- 
ber indices that distinguish quantum Hall states. Quan- 
tum Hall states have non-zero Chern numbers and can 
occur only if time reversal symmetry (TRS) is broken; 
until recently they have been observed only when TRS 
is explicitly broken by an external magnetic field. In 
this article we discuss a class of broken symmetry states, 
first proposed theoretically0-(l] and recently discovered 
experimentally [|1 @], which appear in chirally stacked 
graphene systems and are characterized by spin and val- 
ley dependent spontaneous layer polarization. The aim 
of the present paper is to explain how these states are dis- 
tinguished by their chargejsT-flol]. spin[ll[, and valley [l^ 
quantized Hall conductances, by their orbital magnetiza- 
tions, and by their edge state properties. 

Success in isolating monolayer and few-layer sheets 
from bulk graphite, combined with progress in the epi- 
taxial growth of few-layer samples, has opened up a rich 



new topic[13l Il4| in two-dimensional electron physics. 



Electron-electron interaction effects are most interest- 
ing in ABC-stacked N > 2 layer systems (l5- 17 1. partly 
because [H 0, [l8- 21 j their conduction and valence bands 
are very flat near the neutral system Fermi level. For 
this special stacking order, low-energy electrons are con- 
centrated on top and bottom layers and the low-energy 
physics of a TV-layer system is described approximately 
by a two band model with ±p N dispersion and large as- 
sociated momentum-space Berry curvatures [22j. When 
these band states are described in a pseudospin lan- 
guage, the broken symmetry state is characterized Q by 
a momentum-space vortex with vorticity N and a vortex- 
core which is polarized in the top-or-bottom layers. For 
AB stacked bilayers, for example, interactions lead to 
a broken symmetry ground state [1, 0, 0] with a spon- 
taneous gap in which charge is transferred between top 
and bottom layers. ABC-stacked trilayer graphene has 
even flatter bands and is expected to be even more un- 
stable to interaction driven broken symmetries fl7| . but 
samples that are clean enough to reveal its interaction 
physics have not yet been studied. 

Hall Conductivities and Magnetizations — We discuss the 
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FIG. 1: (Color online) For cases (a-d) the lower panel de- 
scribes the sense of layer polarization for each spin- valley com- 
binations while the upper panel schematically indicates the 
corresponding Hall conductivity contributions, (a) a valley 
Hall insulator with a net layer polarization and a CDW mass 
Acr z ; (b) an anomalous Hall insulator with a valley-dependent 
mass Ar z <r z ; (c) a layer-antiferromagnetic insulator with a 
spin-dependent mass As z <r z ; (d) a quantum spin Hall (or 2D 
Topological) insulator with a valley and spin dependent mass 
term Ar z s z cr z . 



electronic properties of TV-layer ABC-stacked systems in 
terms of ordered state mean-field Hamiltonians of the 
form, 



H N = 



(vop) 



/v 



cos(N(j) f 



■ sin(N<f) p )o-y] + Acr z . (1) 



(-7i) N - 

We have used the notation cos0 p = r z p x /p and sin</) p = 
Py/p where r z = ±1 labels valleys K and K', the two 
inequivalent Brillouin zone corners. The Pauli matri- 
ces in Eq. ([1]) act on a which-layer pseudospin degree-of- 
freedom. In Eq. (fT]) the first term[l{jll7| is the low-energy 
k • p band Hamiltonian for a single valley. Weak remote 
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TABLE I: Summary of spin-valley layer polarizations (T or 
B) and corresponding charge, spin, and valley Hall conduc- 
tivities (e 2 /h units) and insulator types for the three distinct 
states (b-d) with no overall layer polarization on which we fo- 
cus, for a state in which every spin-valley is polarized toward 
the top layer (a), and for a state with partial layer polariza- 
tion. 
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hopping processes have been dropped with the view that 
they do not play an essential role in the broken symme- 
try states |4|. The second term is an interaction- induced 
gap [3, 0, Ho, HH term which defines the direction of layer 
polarization in the momentum space vortex core. For 
each spin and valley, symmetry is broken by choosing a 
sign for A. We have dropped the momentum dependence 
of A because, as we will see, it does not play an essential 
role below. In the mean-field Hamiltonian 2|A| is the size 
of the spontaneous gap, vq is the Fermi velocity in mono- 
layer graphene, and 71 ~ 0.4 eV is the inter-layer hopping 
energy. The p N dispersion is a consequence of the TV-step 
process in which electrons hop between low-energy sites 
in top and bottom layers via high-energy states. 

When spin and valley degrees-of-freedom are taken into 
account, the system has sixteen distinct broken symme- 
try states in which the sign of A is chosen separately for 
(A' t), {K \), (K' t) and (K' |) flavors. We take the view 
that any of these states could be stable, depending on 
details that are beyond current knowledge and might be 
tunable experimentally The sixteen states can be clas- 
sified according to their total layer-polarization which is 
proportional to the sum over spin- valley of the sign of A. 
Six of the sixteen states have no net layer charge transfer 
between top and bottom layers and are likely to be low- 
est in energy in the absence of an external electric field. 
These six states can be separated into three doublets 
which differ only by layer inversion in every spin-valley. 
Thus three essentially distinct states compete for the bro- 
ken symmetry ground state: the anomalous Hall state in 
which the sign of A is valley-dependent but not spin- 
dependent (Acr z — > At z ct z ), the layer-antifcrromagnctic 
state in which A is only spin-dependent (Aer z — > As z er z ) 
and the topological insulator (TI) state in which A is 
both spin and valley dependent (Acr z — >• At z s z ct z ). These 
states are distinguished by their spin and valley depen- 
dent Hall conductivities and orbital magnetizations indi- 
cated schematically in Fig. Q] and summarized in Table Q] 

Berry Curvatures — The three broken symmetry states 
on which we focus are distinguished by the signs of the 



Berry curvatures [22j contributions from near the K and 
K' valleys of f and \. spin bands; we note that the Berry 
curvatures are non-zero only when inversion symmetry is 
spontaneously broken. Using the Berry curvatures, we 
evaluate the orbital magnetizations and Hall conductivi- 
ties of all three states. For momentum-independent mass 
A the Berry curvature of the TV-layer chiral model is 



^4 a) (P>Tz,S z ) 
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dh\\ 
dp 



(2) 



where symbol a = +(— ) denotes the conduction (va- 
lence) band, and the transverse and total pscudospin 

fields are h\\ = (v p) N / 7 f _1 and h t = J X 2 + h\. The 



orbital magnetic moment carried by a Bloch electron|22j 
is m for a two-band model with particle- 

hole symmetry. For the chiral band model 
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where m c is the electron mass and /xb is the Bohr magne- 
ton eh/2m c . Like the Berry curvature the orbital magne- 
tization changes sign when the valley label changes and 
when the sign of the mass term (the sense of layer polar- 
ization) changes, i.e. both are proportional to r z sgn(A). 
The orbital magnetization is however independent of the 
band index a. As illustrated in Fig. [3J in the case of 
|A| = 10 meV, the orbital magnetic moment close to each 
Dirac point has a symmetric sharp peak at which individ- 
ual states carry moments twenty times larger than fiB- 
The state in which A — > Xt z therefore has overall orbital 
magnetization and broken time reversal symmetry, even 
though it docs not have a finite spin-polarization. Inte- 
grating over the valence band, we obtain a total orbital 
magnetization per area ~ (NAm e /2irh 2 ) ln(7i/|A|)/iB, 
that is ~ 0.002/ie per carbon atom for |A| = 10 meV. 

In the presence of an in-plane electric field, an electron 
acquires an anomalous transverse velocity proportional 
to the Berry curvature, giving rise to an intrinsic Hall 
conductivity 0, [22|. Using Eq. ([3]), we find that the in- 
trinsic Hall conductivity contribution from a given valley 
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FIG. 2: (Color online) The magnitude of orbital magnetic mo- 
ments carried by individual states versus in-plane momentum, 
for each spin and valley flavor in ABC graphene iV-layers. 
Here the moments are in units of /ib and |A| = 10 meV. 
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and spin is 

where ht(pF) is the total pseudospin field at the Fermi 
wavevector. Provided that the Fermi level lies in the mass 
gap, each spin and valley contributes Ne 2 /2h to the Hall 
conductivity, with the sign given by T z sgn(A). 

In Fig. QJa) we consider the case in which each spin- 
valley is polarized in the same sense. The total Hall con- 
ductivity is then zero for both spins, with the K and 
K' valleys making Hall conductivity and magnetization 
contributions of opposite sign, preserving time reversal 
symmetry. This phase can be viewed as having a valley 
Hall effect and, even though it does not break time- 
reversal symmetry, we argue later that this designation 
has physical significance. 

As shown in Fig. (Ub), the case Acr z — > Ar z cr z im- 
plies Hall conductivity and orbital magnetization con- 
tributions of the same sign for each spin and valley. This 
state breaks time reversal symmetry but its spin density 
is surprisingly is everywhere zero. The total Hall con- 
ductivity has the quantized value 2Ne 2 /h. Similarly, the 
orbital magnetic moment has the same sign for all flavors. 
We refer to this state as the quantized anomalous Hall 
state. In addition to its anomalous Hall effect, this state 
has a substantial orbital magnetization. The anomalous 
Hall states is probably most simply identified experimen- 
tally by observing a v = 2N QHE which persists to zero 
magnetic field. 

For Acr z — > A,s z er z , depicted in Fig. [He) the two spins 
have valley Hall effects of opposite sign, and the two lay- 
ers have spin-polarizations of opposite sign. This layer- 
antiferromagnetic state has broken time reversal symme- 
try and opposite spin-polarizations on top and bottom 
layers. 

Fig. [Tf d) describes the third type of state with effective 
interaction Acr z — > At z s z ct z . This state does not break 
time reversal invariance and instead has anomalous Hall 
effects of opposite signs in the two spin subspaces, i.e. a 
spin Hall effect. Neither the top nor the bottom layer 
has spin or valley polarization. Quite interestingly if we 
only consider one layer, there are both spin Hall and 
valley Hall effects, however, the orientations of the Hall 
currents in the top and the bottom layers are the same 
for the spin Hall effects but opposite for the valley Hall 
effects. 

Table U includes as well the case in which one flavors 
polarizes in the opposite sense of the other three; charge, 
valley, and spin Hall effects coexist in this state which 
can be favored by a small potential difference between 
the layers. 

Edge States — The physical significance of spontaneous 
charge, valley, and spin anomalous Hall effects is illus- 
trated in Fig. [3] Graphene has very weak spin-orbit in- 
teractions, which in our case we ignore altogether. Fig. [3] 
compares the edge electronic structure of N = 1,2,3 
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FIG. 3: (Color online) Intra-valley and inter-valley edge states 
in chirally-stacked graphene systems. (a)(d) for a single layer, 
(b)(e) for a bilayer and (c)(f) for a trilayer. To visualize 
the edge states, the intralayer and interlayer nearest neigh- 
bor hoppings are chosen as 70 = I and 71 = 0.3, respectively; 
A = 0.25 for Fig.(a,b,c) and A = 0.3\/3t z for Fig.(d,e,f). 



spinlcss models with a quantized anomalous Hall effect 
(i.e. with opposite layer polarizations at two valleys) 
and with a quantized valley Hall effect. The states with 
anomalous Hall effects have N topologically protected 
robust chiral edge states associated with the QHE, as 
shown in Fig. [3Jd)(e)(f). The edge state structure asso- 
ciated with the valley Hall states is more interesting. In 
the N = 1 valley Hall state the Hall conductivity contri- 
bution associated with each valley is 1/2 in e 2 /h units; 
the full unit of Hall conductance requires the two val- 
leys to act in concert. Because they act in opposition 
in the valley Hall state, there is no edge state, as shown 
in Fig. EJa). For TV = 2 on the other hand, each val- 
ley contributes a full quantum Hall effect, and as we see 
in Fig. [3jb) we find two chiral edge states with opposite 
chirality, one associated with each valley. For N = 3 
depicted in Fig. [3jc) , the additional half quantum Hall 
effect from each valley is insufficient to produce a new 
chiral edge state branch. In general we expect [N/2] 
chiral edge state branches at each valley in an W-layer 
stack. Of course valley Hall edge states are topologically 
protected only when the edge-direction projections of K 
and K 1 valleys are not coincident and inter-valley scat- 
tering due to disorder is absent. Nevertheless, we expect 
robust edge states to be present in valley Hall states, as 
found jlH in numerical studies of valley Hall states in- 
duced by an external electric field without interactions. 

Discussion — At the level of continuum-model mean-field 
theory Q, the three charge balanced states we have dis- 
cussed are degenerate. In addition to breaking inver- 
sion symmetry, each breaks two of three additional sym- 
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metries; time reversal (T), spin rotational invariance 
(SU(2)), and the valley Ising symmetry (Z 2 ). The TI 
state preserves only T, the AH phase preserves only spin- 
rotational invariance, and the AF state has Z 2 symmetry. 
Both TI and AF phases break the continuous SU(2) sym- 
metry and therefore Goldstone modes emerge [10J. The 
actual ground state is dependent on subtle correlation 
and microscopic physics issues that are beyond the scope 
of this paper. We note however that it might be possi- 
ble to induce transitions between different possible states 
using external fields. For example, the energy of the 
quantized anomalous Hall state will be lowered by a per- 
pendicular external magnetic field because it has a large 
orbital magnetization. The fully layer polarized state 
will be favored by an external electric field which pro- 
duces a potential difference between the layers. Increas- 
ing the magnetic field further results in quantum Hall 
ferromagnetism 24-2(|. Recent experiments]^, Q in bi- 
layers appear to provide definitive proof that the ground 
state at very weak external magnetic fields is the quan- 
tized anomalous Hall state. 

The quantum spin Hall effect we discuss in this paper 
is in several respects different from that discussed in the 
well known papers [EEI which foreshadowed the identi- 
fication of topological insulators, (i) The quantum spin 
Hall effect is driven by broken symmetries produced by 
electron-electron interactions, rather than by spin-orbit 
interactions ll| which we neglect. The effective spin- 
orbit coupling At z s z <t z due to electron-electron interac- 



tions can be 10 4 times larger than the intrinsic one[27 
(ii) Unlike the previous interaction induced TI phase 
which appears only at finite interaction strengths, here 
the instability to the TI phase is present even for weak 
interactions, (iii) The broken symmetry occurs only for 
N >2, rather than in the single-layer systems [HI (iv) 
Our states are also distinguished topologically, since they 
are characterized by Chern numbers which can have any 
integer value, rather than by a Z 2 label. Of course, only 
N-odd layers are strong TPs, because the helical edge 
modes are likely to localize in a TV-even system due to the 
backscattering process allowed by T28[. Even though 
spin-orbit interactions are very weak in graphene[27|, our 
more elaborate Chern number state classification is there- 
fore not strictly valid once they are included. In this case 
the quantum spin Hall effect will no longer be precisely 
quantized even for odd N. 
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